Mixed Bino-Wino-Higgsino Dark Matter 
in Gauge Messenger Models 



Kyu Jung Bae^, Radovan Dermisek^, Hyung Do Kim^ and 
Ian- Woo Kim^ 

^ School of Physics and Astronomy, Seoul National University, Seoul, Korea, 151-747 
^ School of Natural Sciences, Institute for Advanced Study, Princeton, NJ 08540 

E-mail: hdkim@phya.snu.ac.kr, dermisek@ias.edu, 
iwkimOphy a . snu . ac. kr , baek j 8 1 @phy a . snu . ac. kr 

Abstract. Almost degenerate bino and wino masses at the weak scale is one of unique 
features of gauge messenger models. The lightest neutralino is a mixture of bino, wino 
and higgsino and can produce the correct amount of the dark matter density if it is 
the lightest supersymmetric particle. Furthermore, as a result of squeezed spectrum 
of superpartners which is typical for gauge messenger models, various co-annihilation 
and resonance regions overlap and very often the correct amount of the neutralino 
relic density is generated as an interplay of several processes. This feature makes the 
explanation of the observed amount of the dark matter density much less sensitive to 
fundamental parameters. We calculate the neutralino relic density assuming thermal 
history and present both spin independent and spin dependent cross sections for the 
direct detection. We also discuss phenomenological constraints from b —^ and muon 
f/ — 2 and compare results of gauge messenger models to well known results of the 
mSUGRA scenario. 
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1. Introduction 

Models with weak scale supersymmetry are some of the most attractive candidates 
for extensions of the standard model (SM). Among them the minimal supersymmetric 
standard model (MSSM) is popular due to its minimality and simplicity. Smallness 
of the weak scale compared to the Planck scale is nicely explained by the smallness 
of supersymmetry breaking and the three different gauge couplings meet at the grand 
unification (GUT) scale, 2 x 10^^ GeV, which is close to the Planck scale. Furthermore, 
assuming R-parity, the lightest supersymmetric particle (LSP) is stable and provides a 
reason for the existence of dark matter. 

The lightest neutralino, being weakly interacting, neutral and colorless, and 
appearing as the LSP in a large class of SUSY breaking scenarios including the most 
popular one, mSUGRA, is especially good candidate, since it naturally leads, assuming 
thermal history, to a dark matter density I^dm^^ ~ 1 [H |2]. This observation was 
certainly a major success of supersymmetry. However, the precisely measured value of 
the dark matter density, flDuh"^ ~ 0.1 ± 0.01 [3] together with direct search bounds 
on superpartners tightly constrain supersymmetric models and explaining the correct 
amount of the dark matter density while evading all experimental constraints on 
superpartners is no longer trivial. For example, a bino-like lightest neutralino which 
is typical in the mSUGRA scenario usually annihilates too little which results in too 
much relic density. The bulk region of mSUGRA scenario where neutralino annihilation 
is further enhanced by t-channel exchange of relatively light sleptons and the correct 
amount of dark matter density can be obtained has been highly sqeezed. Indeed, when 
the neutralino is mostly bino and Mi/mg^ < 0.9, the correct relic density constrains 
i^en < 111 GeV at 95% CL [1]. The limits on the Higgs boson mass and b s-j 
independently disfavor the bulk region. In the region with small fi term neutralinos can 
efficiently annihilate via their higgsino components. This region extends along the line 
of no EWSB and is refered to as the focus point or hyperbolic branch region. Remaining 
regions are the regions where special relations between independent parameters occure 
and the neutralino relic density is further reduced by either co-annihilation with other 
superpartners, e.g. the stau co-annihilation region in mSUGRA when stau mass is very 
close to the neutralino mass, or by the CP odd Higgs boson resonance when the mass 
of the CP odd Higgs boson is close to twice the mass of the lightest neutralino. These 
regions require a critical choice of parameters in the sense that the predicted value of 
the dark matter density is highly sensitive to small variations of parameters [1]. 

The lightest neutralino in gauge messenger models is typically mostly bino with a 
sizable mixture of wino and higgsino. The wino and higgsino components enhance the 
anihilation of the lightest neutralino and the correct amount of the dark matter density 
is obtained without relying on critical regions of the parameter space. The virtue of 
the lightest neutralino being a mixture of bino and wino was recognized in studies of 
unconstrained MSSM [H El [3, [H] . As already discussed, the bino-like neutralino typically 
leads to too large relic density. On the other hand both wino-like and higgsino-like LSPs 
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annihilate too fast and the correct amount of is obtained only if they are very heavy, 
m^o ~ 1 TeV for higgsino-like and m^o ~ 2.5 TeV for wino-like neutralino. Obviously 
the lightest neutralino which is a proper mixture of bino, wino and higgsino can lead 
to the correct amount of the dark matter density while avoiding all experimental limits 
and being fairly light. The only problem is that this situation typically does not happen 
in widely studied SUSY breaking scenarios. For example, in models with universal 
gaugino masses at the GUT scale, e.g. mSUGRA, the ratio of bino and wino masses 
at the weak scale is 1:2 and the sizable bino-wino mixing is not possible. However, in 
gauge messenger models the sizable bino-wino mixing is a built-in feature. The bino 
and wino masses are generated with the ratio 5:3 at the GUT scale which translates 
to the ratio 1:1.1 at the EW scale. The bino and wino masses are almost degenerate 
and thus, besides sizable mixing, also the chargino co-annihilation is always present and 
playes an important role. 

The gauge messenger model is independently well motivated [9]. The same field 
which breaks SU{5) to SU{3) x SU{2) x U{1)y, the symmetry of the standard model, is 
also used to break supersymmetry. The heavy X and Y gauge bosons and gauginos 
play the role of messengers of SUSY breaking||| All gaugino, squark and slepton 
masses are given by one parameter and thus the model is very predictive. Besides 
already mentioned non-universal gaugino masses at the GUT scale also the squark and 
slepton masses squared are non-universal and typically negative with squarks being 
more negative than sleptons. This feature leads to squeezed spectrum at the EW 
scale (the heaviest superpartner has a mass only about twice as large as the lightest 
one). Negative stop masses squared at the GUT scale are partially responsible for large 
mixing in the stop sector at the EW scale which maximizes the Higgs mass and reduces 
fine tuning in electroweak symmetry breaking [T2j. Assuming no additional sources 
of SUSY breaking the gravitino is the LSP and then, depending on tan/3, stau or 
sneutrino is the next-to-lightest SUSY particle (NLSP). However, masses of the lightest 
neutralino, sneutrino, stau, and stop are very close to each other and thus considering 
small additional contributions to scalar masses, e.g. from gravity mediation the size of 
which is estimated to be of order 20% - 30% of gauge mediation, neutralino can become 
the LSP. In that case we can utilize the bino-wino-higgsino mixed neutralino feature 
of gauge messenger models to explain the correct amount of the dark matter density. 
Athough there is no necessity to rely on special resonance or co-annihilation scenarios, 
due to the squeezed spectrum of gauge messenger models, these special regions overlap 
and very often the correct amount of the relic density is generated as an interplay of 
several processes. This feature makes obtaining the correct amount of the dark matter 
density much less sensitive to fundamental parameters. 

In this paper we consider the lightest neutralino of gauge messenger models as a 
candidate for the dark matter of the universe. In Sec. [2] we review basic freatures of 
gauge messenger models. We discuss neutralino dark matter in mSUGRA scenario in 

X This idea was suggested in Ref. [10] motivated by Ref . [11] . 
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more detail in Sec. [3] which will be useful when comparing results of gauge messenger 
models. In this section we also outline procedure used to obtain results and summarize 
experimental constraints used in the analysis. Results of neutralino relic density in gauge 
messenger models are presented in Sec. Hltogether with the discussion of constraints from 
b —>■ s'j and muon g — 2. We also give predictions for direct dark matter searches. Finally, 
we conclude in Sec. [51 For convenience, formulae for the composition of the lightest 
neutralino in gauge messenger models which are used in the discussion of results are 
derived in the Appendix. 



2. Gauge Messenger Model 

Let us summarize basic features of gauge messenger models introduced in Ref. [9]. 
The simple gauge messenger model is based on an SU(5) supersymmetric GUT with 
a minimal particle content. It is assumed that an adjoint chiral superfield, S, gets a 
vacuum expectation value (vev) in both its scalar and auxiliary components: (S) = 
(S + X diag(2, 2, 2, — 3, — 3). The vev in the scalar component, S ~ Mg, gives 

supersymmetric masses to X and Y gauge bosons and gauginos and thus breaks SU(5) 
down to the standard model gauge symmetry. The vev in the F component, Fs, splits 
masses of heavy gauge bosons and gauginos and breaks suppersymmetry. The SUSY 
breaking is communicated to MSSM scalars and gauginos through loops involving these 
heavy gauge bosons and gauginos which play the role of messengers (the messenger 
scale is the GUT scale). The gauge messenger model is very economical, all gaugino 
and scalar masses are given by one parameter, 

and it is phenomeno logically viable [9]. 

A unique feature of the gauge messenger model is the non-universality of gaugino 
masses at the GUT scale. The bino, wino and gluino masses are generated with the 
ratio 5:3:2 at the GUT scale: 

Ml = 10 MsusY, (2) 

M2 = 6 MsusY, (3) 

M3 = 4MsusY. (4) 

As a consequence, the weak scale bino, wino and gluino mass ratio is approximately 
1:1.1:2. 

Similarly, soft scalar masses squared are non-universal and typically negative at 
the GUT scale. They are driven to positive values at the weak scale making the model 
phenomenologically viable. Negative stop masses squared are a major advantage with 
respect to the electroweak symmetry breaking which requires less fine tuning and at the 
same time avoids the limit on the Higgs boson mass by generating large mixing in the 
stop sector [12]. In the gauge messenger model the GUT scale boundary conditions for 
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squark and slepton masses of all three generations are given as: 




= -4M|usY, 




(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



m^, = -6M|usY, 
ml = -3M|usY, 
ml. = +6M|usY. 



For completeness, the soft tri-linear couplings are given by: 
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and the same results apply to soft tri-linear couplings of the first two generations. 
The soft SUSY breaking parameters given in Eqs. (I2j) - ( fT3!) correspond to the simple 
SU(5) gauge messenger model with minimal particle content. For soft SUSY breaking 
parameters in extended models see Ref. [9]. 

Gauge mediation does not generate the and 5/i terms and they have to be 
introduced as independent parameters. As we discuss later, they can be generated by 
gravity mediation through Giudice-Masiero mechanism [T3j. The absolute value of /i is 
fixed by requiring EWSB with the correct value of Mz and thus only the sign(/i) can be 
chosen arbitrarily. The other parameter, can be replaced by tan/3 = Vu/vd- Thus 
the simple gauge messenger model has two continuous and one discrete parameters: 



Furthermore, constraints on muon anomalous magnetic moment favor the sign of (i to 
be the same as the sign of the wino mass which in our notation is positive. 

An example of the spectrum of the gauge messenger model is given in Fig. [H For 
comparison we also give a typical spectrum of the mSUGRA scenario in the same figure. 
The mass ratio of the gluino and the lightest neutralino is about 2 in the simple gauge 
messenger model while it is about 6 in the mSUGRA. Assuming no additional sources 
of SUSY breaking the gravitino is the LSP (with the mass of order the EW scale) and 
then, depending on tan /3, stau or sneutrino is the NLSP. However, as is we can see 
in Fig. [H masses of the lightest neutralino, sneutrino, stau, and stop are very close 
to each other and thus considering small additional contributions to soft masses, e.g. 
from gravity mediation or D-term contributions from breaking of f/(l) contained in an 
extended GUT like SO(IO) or E(6), neutralino can become the LSP. In that case we 
can utilize the bino-wino-higgsino mixed neutralino feature of gauge messenger models 
to explain the correct amount of dark matter. 

§ Neglecting mixing in the stau sector, sneutrino would be the NLSP due to the D-tcrm contribution 
which is negative for sneutrino and positive for stau. The mixing in the stau sector is enhanced by 
tan/3 and for tan/3 > 15 stau becomes lighter than sneutrino [9]. 



MsusY, tan/5, sign(/i). 



(14) 
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tan|3 = 10, M =80GeV,n>0 

SUSY 



d,0,0. 



M£2 



(a) Simple Gauge Messenger Model 



tan p = 10, = M,,2 = SOOGeV, Aj= 0, ^l > I 




(b) mSUGRA 



Figure 1. The spectrum of the simple gauge messenger model for tan/3 = 10 and 
AfsusY = 80 GeV (a) and the spectrum of mSUGRA for tan /3 = 10, mo = M1/2 = 800 
GeV, A = (b). The parameters are chosen such that the lightest neutralino mass is 
the same in both cases. 
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Since the messenger scale is the GUT scale, and the gauge mediation is a one loop 
effect, the naively estimated size of gravity mediation induced by non-renormalizable 
operators (suppressed by Mpi) is comparable to the contribution from gauge mediation. 
Gauge mediation contribution is given by Msusy, 

^ (15) 



Mqut 

where C represents group theoretical factors appearing in Eqs. ([2]) - ( |T3l) and the 
contribution from gravity mediation is 

'wl 

For a typical C ~ 5 — 10 and A of order one we find that the gauge messenger contribution 
is about 5 times larger than the contribution from gravity mediation. 

Considering the contribution from gravity, the /i and Bfi terms can be generated [T3] 
together with additional contributions to soft masses of the two Higgs doublets which 
we parameterize by ch^ and c/^^ so that soft masses of the two Higgs doublets at the 
GUT scale are given as: 

mjj^ = - 3M|usY + CH„M|usY, (17) 

mj,^ = - 3M|usY + chMvsy^ (18) 

In addition we also consider a universal contribution to squark and slepton masses which 
we parameterize by cq so that, e.g., 

m| = - llMlusY + coMlusY, (19) 

and similarly for other squark and slepton masses in Eqs. (jSD ^ ©• Thus in the most 
general case the parameter space of gauge messenger models we consider is given by five 
continuous parameters and the sign of the /i term: 

Msusy, tan/3, Cq, ch„, c^^, sign(/i). (20) 

Small contribution from gravity mediation, cq > 5, is enough to make neutralino lighter 
than sneutrino or stau in most of the parameter space. Neutralino is then the LSP or 
NLSP depending on the gravitino mass. Making gravitino heavier is not problematic 
and it can be done assuming other sources of SUSY breaking which do not contribute 
to soft SUSY breaking terms of the MSSM sector. In the next section we consider 
neutralino LSP as a candidate for dark matter. 



3. Neutralino dark matter in mSUGRA 



In the mSUGRA scenario, or in general in any model with universal gaugino masses at 
the GUT scale, the lightest neutralino is a mixture of bino and higgsino. The bino-like 
neutralino typically has a very small annihilation cross section and can not annihilate 
efficiently. As a consequence, if neutralino is the LSP it gives too large relic density and 
thus most of mSUGRA parameter space is ruled out by WMAP data. Representative 
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Figure 2. Slices tlirough tlie parameter space of mSUGRA in mo ~ M1/2 plane for 
tan/3 = 10 (up) and tan/3 = 50 (down) with ^0 = and /i > 0. Blue dots represent 
the region in which the ncutralino relic density is within WMAP range. Shaded regions 
are excluded by various constraints. 
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slices through mSUGRA parameter space are shown in Fig. [2] for tan /3 = 10 and 
tan/? = 50. The white region represents allowed region after various constraints are 
imposed (these are indicated in the plots and will be discussed later) and the blue dots 
represent the region in which the neutralino relic density is calculated to be within the 
WMAP range. In several specific regions of parameter space the neutralino relic density 
is further reduced and these regions are compatible with WMAP data. In the "bulk 
region" (small niQ and M1/2) neutralino annihilation is enhanced by t-channel exchange 
of sleptons. This region is however disfavored by the limit on the Higgs boson mass and 
b s-f. Contrary to the bino-like neutralino, the Higgsino-like neutralino annihilates 
too efficiently and the relic density turns out to be smaller than the WMAP value. When 
the bino mass and the fi term are comparable, sizable mixing is possible and the correct 
relic density is obtained. This is the region in Fig. [2] for larger mo which goes along the 
line where EWSB is no longer possible (or chargino mass limit). In the region where 
the stau mass is close to the neutralino mass the neutralino relic density is reduced 
by co-annihilation with stau. For tan/3 = 10, only tiny stau co-annihilation region is 
available. Bulk region and the focus point region is excluded by direct search bound 
on the lightest Higgs and chargino and the muon anomalous magnetic moment. It is 
indeed the case that the most of the parameter space producing the correct dark matter 
density is already ruled out. For tan/3 = 50, in addition to the stau co-annihilation 
region, the funnel region (pseudoscalar Higgs resonance) appears and also large part of 
mixed Higgsino region (focus point) is not ruled out for larger mo. 

In order to compare the result with gauge messenger models that we will discuss 
later, we fix the ratio of mo/Mi/2 and present a slice through mSUGRA parameter space 
in mo (or M1/2) - tan/? plane in Fig. [3l The region producing the correct relic density 
exists only for large tan /? (tan /? > 45) which is due to the pseudoscalar Higgs resonance. 

3.1. Experimental constraints and procedure used 

Before we discuss results for gauge messenger models let us summarize the procedure 
we use to calculate the neutralino relic density and experimental constraints we employ. 
We obtained our results using SOFTSUSY 2.0.8 [H] for the renormalization group 
evolution of soft SUSY breaking parameters from the GUT scale to the EW scale and 
for calculation of SUSY spectrum. The mass of the lightest CP even Higgs boson is 
calculated with FeynHiggs 2.4.1 [15]. Indirect constraints from b 57, muon g — 2 and 
B — > , and the neutralino relic density are obtained using micrOMEGAs 2.0 [16] 
and the direct detection rates are calculated using DarkSUSY 4.1 [T7] . 
The WMAP result for the dark matter density is [3]: 

^Divi/i^ -0.113 ±0.009. (21) 

In our plots we consider 2cr bounds for the neutralino relic density, 0.09 < Q-Duh'^ < 0.13, 
to be in agreement with the observed dark matter density. This region is represented 
by blue bands in plots. 
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Figure 3. A slice tlirough tlie parameter space of mSUGRA in tan/3 - uiq = M1/2 
plane for = and /.j > 0. Blue dots represent the region in which the neutralino relic 
density is within WMAP range. Shaded regions are excluded by various constraints. 



For B{b — > 57) we consider the allowed range to be 2.3 x 10~^ < B{h 57) < 
4.7 X 10~^ which is obtained by summing the experimental and theoretical error linearly 
and taking the 2a range [HI [T9] . 

Muon anomalous magnetic moment might be the only indirect evidence for the 
presence of new physics at around the weak scale. Recent experimental value of 
Qfj, = {g — 2)^/2 from the Brookhaven "Muon g-2 Experiment" E821 [20] is 

a;fP = (11 659 208 ± 5.8) x 10"^°. (22) 

The standard model prediction contains QED, EW and hadronic parts. As a result 
of undertainties in hadronic contribution, we quote results of two groups for Aa^j = 
a^^P — a'i^ where a^^ stands for the theoretical prediction of the standard model. From 
results of Refs. [21] and [22] we have 

Afl/, = (31.7±9.5) X 10"^°, (23) 

which indicates a 3.3a deviation from the standard model. In order to explain the 
experimental result within 2a, we need a contribution from new physics Aa^j > 
13 X 10~^°. On the other hand, from results of Refs. [23| and [2l|, we have 

Aa^ = (20.2 ± 9.0) x 10^^°, (24) 

which indicates a 2.1cr deviation. In this case we need Aa^^ > 2 x 10~^° if we allow for 
2a variation. Both groups calculated the hadronic contribution using e~^e~ data. The 
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T decay data has not been used because of the uncertainties related to isospin breaking 
effects. By combining these two results [25l[26], we get 

Aa^ = (25.2 ±9.2) X 10"^°, (25) 

which indicates a 2.7a deviation from the standard model. A contribution from new 
physics Aa^ > 7 x 10"^'' is necessary in this case to agree with data. 

In our plots, we draw all three 2a bounds, Aa^ = (2, 7, 13) x 10~^°. As we neglected 
r decay data, we take the most conservative bound, Aa^ = 2 x 10~^°, to constrain the 
parameter space. For illustrative purposes we add two dashed lines corresponding to 
Aa^ = 7 X 10-1° and Aa^ = 13 x 10"^°. 



4. Neutralino dark matter in gauge messenger models 

The squeezed spectrum of gauge messenger models makes the discussion of the 
neutralino relic density very complex. Various regions with correct relic density which 
are usually well separated in scenarios with highly hierarchical spectrum are overlaping 
here and often there is no single process that would be crucial for obtaining the correct 
amount of the neutralino relic density. 

The lightest neutralino in gauge messenger models is typically mostly bino with 
a sizable mixture of wino and higgsino. In order to understand the dependence of 
the neutralino relic density on fundamental parameters it is important to know the 
composition of the lightest neutralino. The formulae for wino and higgsino components 
of the lightest neutralino mass eigenstate are derived in the Appendix and for tan (3 > 10 
they can be writen as: 

Nn ^ 1, (26) 
M| sin 2611^ 
26(/x2 - Ml 
fiMz sin 9w 



^12 - - 77777^ (27) 



'13 



/i2 - Ml 



2 



(28) 



- - - „2_.^2 > (29) 



Ml sin Ow 
/i2 - Ml 

where e is defined as M2 = Mi(l + e). The bino/wino mass ratio is fixed in the gauge 
messenger model. As Mi/gf is RG invariant at the 1-loop level, this ratio at the EW 
scale is 

Mi(Mz) 5 ^ Mi{Mgxjt) 

) - = -tan9w ) r — 0.9, (30) 

M2(Mz) 3 ^M2(Mgut) ^ ^ 

which means e ^ 0.1. From the above equations we see that the wino and higgsino 

mixing is sizable unless the ratio M^/yU is too small. For fi > Mi, the down type Higgs 

component, A^i4, is larger than the up type Higgs component, N13. The bino-wino 

mixing, A'^12 is suppressed compared to the bino-higgsino mixing by Mz/n < 1. This is 

the reason why the bino-wino mixing is negligible in most of SUSY breaking scenarios. 

However, in gauge messenger models the mixing is enhanced by 1/e thanks to near 
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degeneracy of bino and wino. As a result, the lightest neutralino in gauge messenger 
models is mostly bino with sizable and comparable wino and higgsino components. 

Results for the neutralino relic density in gauge messenger models are given in 
Figs, m - [71 We start the discussion with Fig. H] in which we present the results 
for simple gauge messenger model with additional contribution to squark and slepton 
masses, Cq = 10 (up) and Cq = 20 (down). Additional contribution, e.g. from gravity, 
at this level is enough to push masses of all squarks and sleptons above the neutralino 
mass in a large region of the parameter space. Increasing cq shrinks the region of stau 
(N)LSP and opens up the region with neutralino (N)LSP. Blue dots represent the region 
in which the neutralino relic density is within WMAP range. The top part of the blue 
band corresponds to the region of stau co-annihilation. This is easy to understand 
because the blue band stretches along the line dividing the neutralino and stau (N)LSP 
regions. The bottom part of the blue band is due to the CP odd Higgs resonance which 
is not obvious from the plots but it will become clearer in later discussion. In most of the 
region for cq = 10 the co-annihilation with stop is also important and it is the dominant 
process in the region where the two bands meet. However, this does not mean that 
stop co-annihilation and thus the special choice of cq we made is crucial for obtaining 
the correct amount of the neutralino relic density. For Cq = 20 the contribution from 
stop co-annihilation is no longer significant but the shape of the blue band is very 
similar only shifted to the left, to the region of smaller neutralino mass, in which the 
bino-wino-higgsino mixing and the chargino co- annihilation become important. 

For even larger values of cq, see Fig. |5l the effects coming from the exchange of 
or co-annihilation with squarks and sleptons disappear as squarks and sleptons become 
heavy and the band of the correct relic density is independent of cq. The residual small cq 
dependence comes from the fact that increasing cq influences the renormalization group 
evolution of mfj^ in such a way that the size of the /i term increases which consequently 
reduces the mixture of higgsino and wino in the lightest neutralino. As a result, the 
correct value of the neutralino relic density is obtained with slightly lighter neutralino. 

Let us discuss the neutralino annihilation process in detail for one specific point 
from Fig. [5] with Msusy = 42 GeV and Cq = 60. This point is away from the CP odd 
Higgs resonance and the relic density, Q-ouh'^ = 0.11, reflects the composition of the 
lightest neutralino. The lightest neutralino is mostly bino with small mixtures of wino 
and higgsinos: 

A^n = 0.95, A^i2 = -0.22, A^ig = 0.18, Ni^ = -0.09. 

The lightest and the next-to-lightest neuralinos and the light chargino are nearly 
degenerate: 

m^o = 167 GeV, 
m^o = 193 GeV, 

m + = 191 GeV. (31) 

The dominant annihilation channel for this point is XiXi ~^ W^W~ which represents 
31% of the annihilation cross section at the freezeout temperature. It is mediated by the 
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Figure 4. The ncutralino relic density in Msusy - tan (3 plane of the gauge messenger 
model with ch^ = ch^ = 0, /i > and cq = 10 (up) and cq — 20 (down). Blue dots 
represent the region in which the neutralino relic density is within WMAP range. 
Regions with too much or too little relic density arc indicated. For convenience, the 
top axis indicates the mass of the lightest ncutralino. Shaded regions are excluded by 
various constraints. 
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Figure 5. The neutralino relic density in AfsusY ^ co plane of the gauge messenger 
model with ch„ = ch^ =0, /.j > and tan/3 =10. . Blue dots represent the region in 
which the neutralino relic density is within WMAP range. Shaded regions are excluded 
by various constraints. 

t-channel exchange of charginos and thus the wino component of Xi plays an important 
role since the light chargino is mostly the wino. Also important channel is XiXi ~^ 
which contributes 24% indicating that the CP odd Higgs mediated s-channel diagram 
makes a contribution even away from the resonance. This is again a consequence of 
the wino and higgsino mixing (the higgsino-bino-A and higgsino-wino-A interactions 
are crucial). The amplitude for this process scales as A^i4(A'^i2 — tanOwNu) and with 

~ —0.2 we see that the wino component enhances this process by ~ 60%. Finally, 
slepton mediated t-channel diagrams contribute less than 10%. 

Chargino co- annihilation is always present since in the gauge messenger model the 
wino (and thus the lightest chargino) is only about 10% heavier than the bino. The 
chargino co-annihilation for this point contributes about 20% to the annihilation cross 
section at the freezeout temperature. It is mediated mainly by the W boson in the 
s-channel which contributes about 10% and also, to a smaller extent, by the charged 
Higgs in the s-channel. 

In summary, the wino and higgsino mixing and the chargino co-annihilation play an 
important role in obtaining the correct amount of the neutralino relic density in gauge 
messenger models in the region with fairly light superpartners (not ruled out by direct 
searches or the limit on the Higgs boson mass). With this knowledge we can continue 
with the discussion of more typical (and more complex) scenarios when additional 
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co-annihilations and/or resonances further enhance the neutrahno annihilation cross 
section. 

In Fig. [6] (up) we study the dependence of the neutralino relic density on the 
additional contribution to the mass squared for To better understand the behavior 
of the neutralino relic density we also plot the dependence of SUSY spectrum on ch^ for 
fixed value of MgusY in Fig. [6] (down). For > 25 the lightest stop mass is very close 
to the lightest neutralino mass and the stop co- annihilation is dominant. The correct 
amount of the relic density is then obtained in an almost horizontal band at ch^ — 30. 
Going to smaller ch^ the difference between stop and neutralino masses is increasing and 
the co-annihilation with stop is no longer important. The CP odd Higgs resonance takes 
over for ch^ — 20 at MgusY = 60 GeV and somewhat smaller ch^ for larger Msusy • The 
second smaller peak is mainly due to co-anihilation with the lightest chargino through 
the charged Higgs resonance which happens when mH+ — m^o + and to a smaller 
extent due to co-anihilation with the second lightest neutralino through the CP odd 
Higgs resonance which happens when niA — m^o + ^xl' ^i'^'^^ lightest chargino 
and the second lightest neutralino are mostly winos these two resonances happen in the 
same region, ch^ ^ at Msusy = 60, and continue to somewhat smaller ch^ for larger 
Msusy- Finally, decreasing ch^ further takes the lightest neutralino away from stop 
co-annihilation and resonance regions and the blue band of the correct relic density is 
almost vertical in this region. The residual ch^ dependence comes from the fact that ch^ 
changes the size of the /i term which then varies the mixture of higgsino and wino in the 
lightest neutralino. The correct amount of the neutralino relic density in this region is 
obtained entirely due to the wino and higgsino mixing and the chargino co-annihilation 
as discussed in the example above. 

The dependence of the neutralino relic density on the additional contribution to 
the mass squared for H^t is given in Fig. [7] (up). The ch^ controls masses of the heavy 
CP even, charged and CP odd Higgs bosons and only negligibly affect everything else. 
Thus the region of the correct relic density is a vertical band except for the CP odd 
and charged Higgs resonances. Finally, in Fig. [7] (down) we chose such values of ch^ 
and ch^ that the CP odd Higgs resonance does not appear. This plot is similar to those 
in Fig. m but now the stau co-annihilation region turns into a vertical band signaling 
independence on tan/3. Similar vertical band appears also in the mSUGRA scenario, 
see Fig. |3l but it is in the region ruled out by direct searches for SUSY and the Higgs 
boson. 

4-1. Discussion ofb^s'-f and muon g — 2 

From Figs. |l]-[7]we see that the limits on B{b 57) are typically the most constraining 
out of all direct and indirect limits. The charged Higgs contribution is additive to the 
standard model contribution and scales as 



B{b — > S7)i^± oc 




(32) 
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c„ = 0, c = 10, tan B = 10, u > 0| 

1Q0 2p0 3Q0 
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Figure 6. Left: the neutralino relic density in Msusy - ch^ plane of the gauge 
messenger model with cq = 10, ch^ = 0, /U > and tan/3 = 10. Blue dots represent 
the region in which the neutralino relic density is within WMAP range. Shaded regions 
are excluded by various constraints. Right: the dependence of various superpartner 
masses on ch^ for the choice of parameters corresponding to the plot on the left with 
Msusy = 60 GeV. 
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MsusY(GeV) 

Figure 7. Left: the neutralino relic density in Msusy - ch^ plane of the gauge 
messenger model with cq = 10, ch^ = 0, /x > and tan/3 = 10. Right: the neutralino 
relic density in Msusy - tan/3 plane of the gauge messenger model with ch^ = —50, 
ch^ = —50, /i > and tan/3 = 10. Blue dots represent the region in which the 
neutralino relic density is within WMAP range. Shaded regions are excluded by various 
constraints. 
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while the chargino-stop loop contributes as 

B{b s-f)i oc . (33) 

The chargino-stop loop contributes with opposite sign compared to the charged Higgs 
diagram if fiAf is negative. In gauge messenger models the charged Higgs is typically 
heavier than stop and the chargino-stop loop dominates the new physics contribution. 
As a result, the predicted branching ratio becomes lower than the Standard Model result 
and the lower bound on B{b — *• sj) plays an important role. 

In the limit when Mi, M2, and m^;^ are approximately equal, which is the case 
in gauge messenger models, and > Mi, M2, the expression of the supersymmetric 
contribution to muon anomalous magnetic moment [27] simplifies to: 

- 327r2 M2(/x2 _ M2) ^'^^^ 

where M represents sneutrino, smuon, chargino or neutralino masses. It can be rewritten 

as 

A<"--13(»V '^V-/'-0-'° (35) 



As a result, we obtain a relation between M and tan/5. In most of the parameter 
space /i is just about twice as large as the lightest neutralino mass and thus we can set 
J^M— ~ - in which case we get 



-M2 — 3 



26 /100GeV\ 



2 



Aa^ X 10^" - - j tan/3. (36) 

Assuming conservative bounds 2 x 10^^° < Aa^ < 50 x 10~^° a discussed in Sec. 13.11 
we can derive the lower and upper bounds on M as a function of tan/?: 

Mower ~ 40Vtan/5 GeV, (37) 

and 



Muppcr ~ 200Vtan/5 GeV. (38) 

For tan /? = 10 we find 130 GeV ^ M < 630 GeV and similarly for tan /3 = 50 we have 
280 GeV ^ M < 1400 GeV. In Figs. |4]-[7]the value of M approximately corresponds to 
the neutralino mass represented by the top axis. It is interesting to note that the indirect 
bound from the upper limit on the muon anomalous magnetic moment is already well 
above the direct search limits on superpartners. 

As a result of the squeezed spectrum of gauge messenger models the limits on 
B{b s'-f) and the muon g — 2 are almost parallel to each other, see Figs. |1]-[71 This 
is a consequence of the SUSY contribution to both processes scaling approximately as 
tan/3/M^. The limits on B{b s'-f) constrain the SUSY spectrum from below while 
the limits on g — 2 constrain the parameter space from above. The allowed parameter 
space is then only a strip in between these two bounds. This is a characteristic feature 
of models with squeezed spectrum. If the required value of Aa^ turns out to be close to 
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the upper range of current estimates most of the parameter space of gauge messenger 
models we considered will be ruled out with only tiny regions remaining, see Figs. |1]-[71 
Interestingly, it is still possible to obtain the correct amount of the dark matter density 
in these tiny regions, see Figs. [6] and [7] (down). 



4-2. Direct Detection of Neutralino Dark Matter 



In this section we calculate the spin dependent and spin independent nuetralino-nucleon 
cross sections in gauge messenger models. Spin independent neutralino-nucleon cross 
section is dominated by (light and heavy) Higgs mediated t-channel diagrams which are 
controled by the higgsino component of the lightest neutralino: 



2/2 4 



14 



m 



H 



(39) 



where = /t, + ^/tg and X^ = fr^ + ^/tg H- 

Substituting A''i3 and Nu from Eqs. (l28l) and (l29i) we find 



a 



/4 4 
9 T^N 



2\2 



Xd tan P 



m 



H 



X, 



nit 



(40) 



S quark exchange diagrams are negligible due to the hypercharge as long as the squark 
masses are comparable to the heavy Higgs mass. Inserting the numbers X^ = 0.144 
and Xd = 0.18 [2], we get the direct detection rate close to the one we obtained using 
DarkSUSY. 

The detection cross sections for points with the correct neutralino relic density from 
Fig. m (up), the gauge messenger model with ch^ = ch^ = 0, /i > and Cq = 10, that 
satisfy all direct and indirect constraints are given in Fig. |H1 In gauge messenger models 
with no additional contribution to Higgs soft masses and only small contribution to 
other scalar masses enough to make them heavier than the lightest neutralino the direct 
detection cross section scales as a^N oc tan^ /3/MgugY for tan/3 > 10. This behavior is 
clearly visible in Fig. [8l The thickness of the line is determined by the allowed region 
for tan/5, in this case 5 < tan/5 < 25, see Fig. HI The predicted cross sections are 
not within the reach of CDMSII [29]. Assuming additional contributions to Higgs soft 
masses allows for a wider range of the higgsino and wino mixing and the range of the 
predicted detection cross sections spreads as is shown Fig. [91 Part of the parameter 
space is within the reach of CDMSII and the whole parameter space of gauge messenger 
models can be explored at Super-CDMS [29] . 



5. Conclusions 



The lightest neutralino in gauge messenger models is mostly the bino with a sizable 
mixture of the wino and higgsino. The wino and higgsino components enhance the 
neutralino annihilation cross section. Furthermore, the splitting between the bino and 
wino masses is at the level of 10% and thus the co-annihilation with the chargino 
is contributing in the whole region of the parameter space. These two features, the 
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Figure 8. Spin dependent and spin independent neutralino-nuclcon cross sections for 
points with the correct neutralino relic density from Fig. |4] (left), the gauge messenger 
model with ch^, = ch^ = 0, /i > and cq = 10, that satisfy all direct and indirect 
constraints. The lines represent the current CDMS limits |28j and expected limits from 
CDMSII [53] for spin independent cross section. 



lightest neutralino being a mixture of the bino, wino and higgsino, and the chargino co- 
annihilation, are sufficient for obtaining the correct neutralino relic density to explain 
WMAP results with fairly light neutralino (and other superpartners) while satisfying all 
the constraints from direct searches for superpartners and the limit on the Higgs boson 
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Figure 9. Spin dependent and spin independent neutralino-nucleon cross sections 
for points with the correct ncutralino relic density satisfying all direct and indirect 
constraints obtained in an extended scan over whole parameter space of gauge 
messenger models discussed in this paper. The lines represent the current CDMS 
limits and expected limits from CDMSII for spin independent cross section. 



mass. 

This is in contrast with scenarios with the usual hierarchical spectrum, e.g. 
mSUGRA, in which the properties of the lightest neutralino (being bino-like) typically 
lead to the correct neutralino relic density in the region which is already ruled out 
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by direct SUSY and Higgs searches or disfavored hj b S7. In mSUGRA-like 
models obtaining the correct amount of the neutrahno rehc density rehes on special 
co-annihilation and resonance regions which are critically sensitive to small variations 
of independent parameters. Due to a large hierarchy in the spectrum these surviving 
strips are typically well separated by large regions of the parameter space ruled out by 
WMAP data. 

In gauge messenger models, as a result of the squeezed spectrum of superpartners, 
various co-annihilation and resonance regions overlap and very often the correct amount 
of the neutralino relic density is generated as an interplay of several processes. For 
example the stop co-annihilation contributes significantly in a large region of the 
parameter space. This can be easily understood from the fact that both stop and 
neutralino masses are mainly controled by the same parameter and as it happens the 
neutralino and stop masses are very close to each other. Varying contributions to scalar 
masses from other sources is only slowly changing this relation. Furthermore, even if 
we increase stop masses by assuming an independent additional contribution, which 
effectively shuts down the stop co-annihilation, the band of the correct neutralino relic 
density only moves to the region with somewhat lighter neutralino which still satisfies 
the limits from direct SUSY and Higgs searches. This feature makes the explanation 
of the observed amount of the dark matter density much less sensitive to fundamental 
parameters. 

In gauge messenger models with no additional contribution to Higgs soft masses 
and only small contribution to other scalar masses enough to make them heavier than 
the lightest neutralino the direct detection cross section is predicted to be in the range 
10~^^ cm^ - 10~^^ cm^ which is not within the reach of CDMSII but can be explored at 
Super-CDMS. 

Some of the results concerning the neutralino relic density in gauge messenger 
models, namely the presence of various co-annihilation regions, originate from the 
sqeezed SUSY spectrum. Therefore we expect similar results for other models derived 
in different contexts which lead to squeezed spectrum, e.g. deflected anomaly mediation 
[30] [31] [32] and mirage mediation [33] [M] [35] [36] • However, the special features of 
the gauge messenger model related to the bino-wino-higgsino mixed dark matter and 
with that associated chargino co-annihilation depend on details of a model and is not 
automatically guaranteed by the squeezeness. 

In conclusion, let us note that both natural EWSB and natural explanation of 
the correct amount of the dark matter density independently disfavor models with 
hierarchical spectrum. Models with squeezed spectrum seem to be favored and thus 
it is desirable to explore their phenomenological and collider predictions. 
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Appendix 



M, 



N 



(41) 



In this appendix we set conventions and derive approximate formulae for the composition 
of the lightest neutralino in the gauge messenger model which are useful for the 
discussion of neutralino relic density. 

The neutralino mass matrix in the basis (B, W, hd, h^) is given by: 

( -Mzci^se^ Mzspse^ \ 

Ms Mzcpce^- -MzSpCe^ 

-MzCfsSg^ MzCfsCg^ -/i 

y MzSfiSg^ -MzSpCe^ -/i 

where s^^- = smOw, cq^ = cos9]y with 9^ being the Weinberg angle (weak mixing 
angle) and similarly S/s = sin/3, C/3 = cos (3 where tan/3 = 

In the gauge messenger model, bino and wino masses are comparable. Thus it is 
convenient to express the wino mass in terms of the bino mass and a small parameter 
describing the difference, 

M2 = Mi(l + e). (42) 

Numerically e ~ 0.09 and it is almost independent of tan/3. Thus, for Mi < the 
lightest neutralino is mostly bino and the splitting between the bino and the wino is at 
the level of 10 %. 

The neutralino mass matrix can be brought to a diagonal form by an orthogonal 
t r ansf ormat ion , 

Mdiag = A^M^viV^, (43) 

where Nij, j = 1,2,3 and 4 represent the mixture of B, W, h^, hu in the lightest 
neutralino mass eigenstate. 

In order to calculate Nij, it is convenient to rotate the neutralino mass matrix to 
a basis in which the lower right 2x2 block is diagonal, 

(Ml - 

M2 

■^MzS0^{si3 + cp) -^Mzce^isp + cp) 

\ -^^zsewisp - cp) -^Mzce^isp-cp) 

which is obtained by an orthogonal transformation, 

M = UMnU^, 



M 



-^Mzsey^{sp + cp) 
-^MzC0^{sp + cp) 





-^Mzse^{sp- cp) \ 

-^MzCB^{sp - Cfj) 





-At 



(44) 
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with 



U 



( 1 










1 

1 






1 

1 



(45) 



The matrix M can be diagonahzed by an orthogonal transformation, 



VMV^. 



(46) 



The advantage of M is that we can treat off-diagonal elements as perturbations and 
calculate eigenvectors (elements of V) using matrix perturbation formalism. Then, the 
mixing matrix (the diagonalization matrix in the original basis) is simply given as 

N = VU. (47) 

In the leading order, neglecting the off-diagonal elements of M, the diagonalization 
matrix V is an identity matrix. When the mass differences between eigenvalues are not 
extremely small, (M2 — Mi)/i^ > MiM|, or equivalently e/i^ > M|, non-degenerate 
perturbation formalism can be applied. At the first order of perturbation theory we 
have: 



where 



Similarly, the second order corrections are given as: 



(2) 



nl 



MlmMrr 
AnlAfit 



For Ml < we find: 

^12 



0, 



(1) 



13 



(1) 



14 



+ 



M31 _ Mz sin 9w (sin /3 + cos P) 
V2ifi - Ml) ' 
Mz sin 6*14/ (sin (3 — cos (3) 



fi- Ml 
M41 



y2(Mi + /i) 



(48) 
(49) 

(50) 
(51) 

(52) 
(53) 



and thus the higgsino component in the lightest neutralino appears at the first order. 
Since v}^^ = it is necessary to calculate the contribution from the next order. This 



contribution is small in general but can significantly alter the result when Mi 
The second order correction is 

M23M31 M24M41 



Mo. 



y{2) 



12 



(M2 - Ml) (Ml - /i) (M2 - Ml) (Ml + /i) ' 



M| sin 26'H/(sin p + cos Pf 



M| sin 26 



w 



[sm P — cos P)" 



4eMi(/i-Mi) ' 4eMi(/i + Mi) ' 

and, since e ~ 0.1, it is comparable to the first order corrections coming from the higgsino 
mass. Therefore, we have a sizable bino-wino mixing in addition to bino-higgsino mixing. 
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The diagonalization matrix V is then approximately given as ~ 1 + V^^^ + V^'^\ 
Finally, we can find the components of the mixing matrix in the original interaction 
basis. Using Eqs. (H7|) and ( 145|) we get: 

Nn ^ 1, (55) 
Nu^vS\ (56) 

11 Mz sin Owifi sin (3 + cos (3) 
^13 = + ^V,, + -^V,, ^ , (57) 

1 Mzsin6'H/(Misin/? + /icos/3) 



'14— 7^*^13 "I 7^ ^14 



(58) 



V2 V2 /i^-Mf 
For tan/5 > 10 these formulae can be further simplified: 

iVii ^ 1, (59) 
M| sin 2ew 



'12 



2e(/i2 - Mf) ' 
jjiMz sin 9w 

/i2 - Ml 
Ml sin Ow 



(60) 



'13 ^ + —2 (61) 



'14 



/i2 - Ml 



(62) 
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